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CERTAIN INTEGRATION FORMULAE USEFUL IN 
NUMERICAL COMPUTATION. 



By S. A. COREY, Hiteman, Iowa. 



In developing a function of a variable by Taylor's, Maclaurin's, or the 
Euler-Maclaurin formulae, the evaluation of the successive derivatives is a 
fundamental and in most cases indispensable operation. In many instances 
the higher derivatives become exceedingly complex and their evaluation cor- 
respondingly difficult. In such cases frequently the labor involved in finding 
results of the required degree of accuracy by these formulae becomes pro- 
hibitively great, and other methods are resorted to, involving usually some 
of the formulae of the Calculus of Finite Differences. The formulae given 
below, some of which are probably new and others old except as to the re- 
mainder term, will, it is believed, prove useful in many such cases, for the 
reason that the use of the higher derivatives is not involved except in the 
remainder term, in which term a roughly approximate value usually suffices. 
In certain cases it will be found desirable to develop the given function to a 
limited number of terms by Taylor's or Maclaurin's formula, to differentiate 
this series the number of times required to get an expression for the deriv- 
ative involved in the remainder term, and then from this expression to de- 
termine roughly the value of the remainder term. Another method of indi- 
cating the degree of accuracy assured in cases where the value of the 
remainder term is not easily determined is to develop the function by two of 
the formulae given below, and then to compare the results obtained. It is 
usually safe to assume that the difference in the results obtained may be 
substituted for the remainder term in either of the two formulae employed. 
It will be noted that the remainder term, as given below, is expressed, in a 
form similar to that frequently employed in Taylor's formula, and that it is 
therefore always possible to compare the accuracy assured with the accuracy 
assured by the use of Taylor's formula. 

From the manner of deriving these formulae it will appear that it is 
not practicable to give all the formulae of like form, obtainable from the 
general formula (c), but the list here given is, probably, sufficiently complete 
to satisfy all ordinary needs. Should others be desired they may, of course, 
be obtained from the general formula. In fact, were it not for the consid- 
erable amount of labor involved in deriving these formulae one at a time, all 
these formulae could be thus obtained and there would be no occasion for 
appending the following rather formidable list. 

These formulae are all based on the so-called Euler-Maclaurin formu- 
la, a semi-convergent formula which has been recently studied by E. Borel, 
E. W. Barnes, * and other European mathematicians, but to which but scanty 

* Borel, Lecons sur les series divergentes. Chap. III. 
Borel, Acta Mathem&tica, Vol. XX, page 360. 

Barnes, Quarterly Journal of Mathematics, Vol. XXXV, pages 175-188. 
Barnes, Proc. London Mathematical Society, Series 2, Vol. 3, page 253. 
Cantor, Geschichte der Mathematik, Vol. 3, pages 635, 661. 
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reference is made in our American mathematical literature. It would seem 
a greater use of this formula in our American text books would be advisable. 

William Fleetfoot Sheppard, in his article on "Mensuration" in the 
11th edition of the Encyclopedia Brittannica, volume XVIII, pages 142-3-4, 
paragraphs 68 to 79 inclusive, gives a very satisfactory reference to this 
class of formulae, to which the reader is referred. In paragraph 71 he gives 
a general method of constructing some of these formulae, but it should be 
carefully noted that the subscripts of C which he uses are determined by a 
system of notation somewhat different from that herein used. In paragraph 
79 he gives another general form involving an expression for the remainder, 
but the use of the formulae in general numerical work is not emphasized, 
and the reader is likely to infer that their use is confined chiefly to the * sub- 
ject under discussion. 

The Euler-Maclaurin formula may be written, 



f(a+x)=f(a) + -^{f (a) +f (a+x)+2[f (a+^)+f (a+ |?) 



+f(a+^)....+/(af^)]}-|3 r |,[/"(a+ a; )-r(<x)] 



•f^L/Ma+x) -/*(«■)].... 



+ ( ~ 1)?t mM2n) ! [/»(«+«> ~/ 2 "(a)] +S-. (a). 



B n being Bernoulli's nth number. 

Expressed briefly and for clearness in what follows, this becomes: 



f(a+x)=F(x, m)+— H— 1+ h+ -?.... + %+S....(b), 

J m m 4 m e m s m 2n 



in which k 2 , k it k 6 , ...., k 2n , are independent of m, and Sis an expression 
for the remainder after the general term. 

If (r+1) values of m be taken, r of the quantities k, say hi, 
&2(i+n, k2u+z), ...., fou+r-i) , may be eliminated between the (r+1) resulting 
equations, and the result expressed in determinant form as follows: 
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f{a+x) = 





M,+Si 

M 2 +S 2 


m ~2i 
m -2i 


m -2«+l) 
m -2(t+l) 


m -2(»+2) 


.... mf 2 «+'^ 1 > 
.... m^ 2K+r - 1) 




M(r+1) +S(r+1) m^ 1} 


wj— 2K+1) 


m — 2W+2) 
'"■(r+l) 


.... '»( r+ i) 




1 

1 


m -2i 

m^ 2i 


m -2W+l) 

W -2U+1) 


m -2(i+2) 
m -2(i+2) 


.... mf 2 "^-^ 
.... m- 2(i+r - 1 




1 


m Tr+l) 


w — 2(t+l) 
"V+D 


otj— 2(i+2) 
m (r+l) 


<*M — 2(i+r— 1) 
.... '^( r+ i) 



• (c), 



where .Mj, M 2 , M(h-d, represent the approximate values of f(a+x), 

obtained by (b) for m x , m 2 , ...., mw+i), respectively, by employing all the 

terms preceding the term —£-, and where S u <S 2 , ..... S^+d, represent the 

corresponding remainder terms, i. e. the remainder after the term 

m 2«+r~i) m eac ^ of the (r+1) equations. 

Splitting up the determinant of the numerator into two other similar 
determinants, one containing the M elements and the other the S elements, 
we have, by letting D represent the determinant, 



f(a+x)=D M /D+D s /D.... (d). 



To rid the indices of the elements of these determinants of the minus 
sign, divide each element in a row by the last element in that row. Now 
reduce the order by taking the co-factors of all the elements in the first 
column. These co-factors have evidently as factors (m 2 u — in i v ), u and v 
taking any value assumed by the subscripts of m, except that of the ele- 
ment, the co-factor of which is involved, and except u=v. Their evaluation 
follows readily. 

The remainders, S, used in determining the remainders in the follow- 
ing formulae, were obtained from the solution of problem 237, Calculus, given 
in the June-July, 1907, number of The American Mathematical Monthly, 
to which reference must be had for an understanding of the method used in 
their derivation. Each remainder term in the following is dependent for its 
value on but one higher derivative, and for that reason the limits of its 
value are fixed, while the form given by Mr. Sheppard is a function of all 
derivatives higher than a certain derivative, and usually becomes divergent 
if a sufficiently large number of terms be taken. 

Instead of retaining in the following formulae the M\, M 2 , 

M( r +i) , used in (c), it is convenient to employ a slightly different notation, 
viz., C m represents the chordal area of the curve y=f (a+x), where the 
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value of y is taken at each of the (to+1) equidistant intervals, 0, x/m, 2x/m, 
Sx/m, ...., (to— l)x/m and x, thus 



C=j-[ /,(o)+ ^ (a+a?) +/ (« +5-) +/' (a+ ?r) +••••+/' (<*+ %r^)]- 



to 



to 



m 



m 



In the formulae involving y, the first extreme ordinate is invariably 
y , and the last extreme ordinate y n , where n has the highest value of the 
subscript given in that formula. The distances of the intermediate ordi- 
nates y r (where o<r<n) and the last ordinate y n from the first extreme or- 
dinate are proportional to their subscripts. The ordinates are therefore not 
all equidistant in cases where the smaller values of m are not all factors of the 
largest value. 

Where C is given more than one subscript, thus C m . t, it is understood 
that to the expression for the chordal area Cm has been added the sum of the 

ki ' -| + etc., in (b), where 



first t terms of the series 



2 ' ™,l' 



TO 



TO* 



m" 



k 2n =+(-l) n ^{[f 2 Ha+x)-P"(a)l 



B n being Bernoulli's »th number. 



f(a+x)dx =j ydx, and R*=j^\.f a "(a+o 1 x)-f**(a+e t x)], 
where both o % and <? 3 have some value greater then zero and less than unity. 



ii=l(4C,-C,)±A.R4 ..... 

4=4(90,-0 ±Ai?4 

4=jV(64C 4 -20C 2 .+O±T 3 <ri?6 • • • 

A=rw(243C 8 — 128C,+5d) ±-frR 9 
i i= ¥1 ^(&92C4-65ttC 8 +896C,.-7O±.0229i2 8 
A=riT(612C 1 ,-175C a +14C 11 -O ±.000,14# 8 
A=A(66C, a -28C e +8C 4 -0 ±.000,099fl 8 . 
A=A(15C e -6C,+C,)±.006,62Be . 
A=rfB (567C 6 — 512C 4 +162C 3 -7C* ) ± .002,55i2 8 
A=idnr(10,368Ce -8,192C 4 +2,025C 3 -C, ) ±.001,93# 8 



(1). 
(2). 
(3). 
(4). 
(5). 
(6). 
(7). 
(8). 
(9). 
(10). 



122 

A=^(6,661C e -3,584C 4 +175C«-2C I )±.003 f 77iJ 8 . . . (11). 

A= B i T (l,296C 6 -567C 8 +112C,-C 1 )±.008,39ie 8 .... (12). 

A=nW34992Ce -32,768C 4 +10,935C 3 -560C,+d) ±.002,74# 10 (13). 

A= I ^ r (4,096C 8 -l I 844C 4 +84C 2 -C 1 )±.002,02B 8 . . . (14). 
A= T? „i OTff (33 ) 554,432C 8 -31,177,872C 8 +7,569, 408C 4 -1,240,029C 3 

+25,872C 2 — llC,)±.000,225i2 12 (15). 

A=rTnmnnr (2,965,984(7, a -l,667,952C e +585,728C 4 -104,247C 8 

+2,288C 2 -C,)±.000,095,6.R 12 (16). 

^ = iirW(l,728C, 2 -945C e +320C 4 -54C 3 +C 2 )±.000,071,2B 10 . (17). 

A= IBT (384C 6 -105C,+C 1 )±.006,05Be (18). 

4=C 6 +M(C 6 -C 3 )-^(C 3 -C,)±.006,05fl 6 .... (18a). 

A=^(248C e -85C,+2C,)±.017,liJe (19). 

A=C e +n(C e -C i )-rls(C 2 -C I )±.0n,lR e .... (19a). 

4= 3 V(64C 8 -20C 4 +C 3 )±.001,34i?e (20). 

^=7WTnnr(l,048,576C, e -348,160C 8 +22,848C 4 -340C 8 + C,) 

±.000,030,2^ (21). 

^=^™<nnnny(110,075,314,176C 2 4 ^64,900,362,240C, , +25,081,937,920C 8 

— 5,178,990,960C e + 202, 076, 160C 4 -9,808, 695C 3 +38,640C 2 — C,) 

±.000,000,134J? 16 (22). 

A= 2 ^(3,200C,o-825C 5 + C 1 )±.000,624B e (23). 

^=^tW^7^(3,648,810,176C, 8 -1,674, 098, 987C 9 +36,837,504C e 

-575,586C 3 +20,672C 2 -11C,) ±.000,000, 736i? 12 . . . (24). 



Formulae involving dy/dx. 
A- T V(16C 2 ,i-Ci.i)±.115i2 e 
A=C 2 ,i+TV(Cu-Ci,i)±.115i? 6 . . 

A=^(81C 3 .i-Ci,i)±.019,li2 e . 
A=C3.i+ ¥ V(C3.i-Ci,i) ±.019,LR 6 
A= 7 k(l,024C 4 ,i-80C 2 ,i+Ci,i)±.005,9i2 8 
^=C4,i+^A(C4.i-C 2 ,i)-^f(C2,i-C 1 ,i) ±.005,9i2 8 - 
A - T ^D (2, 187C 3 ,i - 512C 2 .i+5Ci,i) ± . 019, 6i2 8 . 
A =rshsjs (131, 072C 4 ,i-59, 049C 3 ,i+3, 584C 2 ,i - 7C U ) ± 



002,88i2 1o 



"8877F 



(36,864Ci 2 ,i -3,150C 6 ,i+63C 3 ;i-2C 2 .i) ±.000,002,76# lo 



(25). 
(25a). 

(26). 
(26a). 

(27). 
(27a). 

(28). 

(29). 

(30). 
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A=rir(270C M — 27C M +2Cb)±.000,3892? 8 (31). 

A= in ^r(10 > 206C M -4096C w +729Cs,i-14C u )±.000,15LB 1 . . . (32). 

ii= ¥in jW(873 l 248C w — 1S1 I 072G U +18,225C U — C u ) ±.000,11&R 10 . (33). 

A=-Bvhrs(l18,mC«j.— 28,672C 4 ,i+350C2.i— C u ) ±.000,2572 1 . . (34). 

A=^^(46,656C6,i-5,103C3,i+448C2,i-Cu)±.000,633i2 lo . . (35). 
A= ¥¥T V^(l ) 259,712C 6 ,i-524,288C4,i+98 ) 415C3. 1 -2,240C2.i+Cu) 

±.000,20272 12 (36). 

A=Wt>7s(262,144C8,i-21,504C4.i+336C 2 ,i-Ci.i) ±.000,089,872 10 . (37). 
A= TT¥rT V^TO(2 ) 147,483,648C 8 ,i-l,122,403,392C 6 ,i+121,H0,528C4,i 

-ll,160,261Cj.i+103,488C w -llC u ) ±.000,010,7/2,4 . . (38). 
A=^ ¥ ^ow(429,981 ) 696Ci2 > i-60 ) 046,272C 6 , 1 +9,371,648C 4 ,i 

-938,223G,i+9,152C 2 .i-Ci,i)±.000,002 ) 8i? 14 . . . (39). 

A=Wm (124, 416Ci2.i - 17,'OlOGu +2, 560<7 4 ,i-243C 3 ,i + 2Gu) 

±.000, 001, 6672, 2 (40). 

A=^A^(320,OOOCio.i-20,625C 5 ,i+Ci,i)±.000,014,3i2 8 . . . (41). 



Formulae involving -Jr., -^ and - ' f - 



dx' dx z 

A=tSr(64C w -Cu2)±.078,4B 8 . 

A=C2,2+MC^-Ci, 2 ) ±.078,472 8 

A=7^r(729C M -Ci.2)±.006,03B 8 

A=C 8i 2+7i»(C M -Ci.2)±.006,0aB 8 . 

A=tti W (16384C4,2-320C 2> 2+ Ci, 2 ) ± . 001,3172 x 

A= T7 i ¥ o(19 ) 683Cs,2-2 ) 048C2,2+5Ci,2)±.007,05B 10 

A=C3,2+tWo (^3,2 ^2,2) — TbV^Gw - C1.2) ±.007,0572, 



dx 6 ' 



(42). 
(42a). 

(43). 
(43a). 

(44). 

(45). 
(45a). 



=T5HkTxr(2,097,152C4,2-531,441C3,2+14,336C 2> 2-7Ci,2) ±.000,66972 ! , (46). 

^x^T5(5,308,416Ci2.2-113,400C6,2+567C3.2-8C2.2) 

±.000,000,086,372! , (47). 

A-Ww(9,720C 6 ,2-243C3,2+8C2.2) ±.000,037,972, „ . . . . (48). 
A=C M +Tf*7(C M -C M )- T A T (C M -C\2)±.000 f 0S7,91? 1 o . . (48a). 



A= 



JW (45, 349, 632C 6> 2— 8, 388, 608C 4 , 2 + 885, 735C 3>2 - 8, 960C 2 , 2 



stsstsc 

+Ci, 2 )±.000,027,372 14 (49). 

A=TTratTHTr(l,679,616C M -45,927C M +l,792Cb-C u ) ±.000,07972, , (50). 
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A=TTm^n r (16 f 777,216C M -8H064C<2+l f 344Ci i 2-Ci.2) 

±.000,006,39i2 12 (51). 

^=^wTk7nyo (61,917,364, 224Ci2,2-2,161,665,792C 6 ,2+149,946,368C4.2 



-8, 444, 007C 3 .2+36, 608(72,2 - Ci, 2 ) ±.000,000, 108B x « 
A=^(256C 2 ,3-Ci,3) ±.073,172,0 ..... 
A=C M + I i T (C M -Ci.»)±.073,lfl 1 , .... 
A= T1 5W(6,561C8,8-Ci.8)±.002,53B 1 o 
A=C w +Wmf(Csjt-Cij)±.Q02,B8R l0 
A=nnAnnr(26^144C u -l,280Cta+Cu8) ±.000,39LR, , 
A=C 4>3 +^oV^ [l,279(C 4 .3-C 2 ,3)-(C2.3-Ci,3)] ±.000,39LR, . 
A=TnnnAnnn7 (60, 466, 176C 6 , 3 -413, 343C 3 .3+7, 168C 2>3 - C 1>3 ) 

±.000,013,372,4 



Formulae expressed in terms of ordinates. 

or 
A^-jrliyo+y^+iy,] ±^ s R i (Simpson's first formula) 

A=-g-[(2/ +2/ 3 ) +3(2/i +2/ J] ±*VR 4 (Simpson's second formula) . 

A=|q[7 (2/0 +2/4) +32(2/, +2/ 3 )+122/ 2 ]± 7 ye 6 

A=~[ll(2/o+2/6)+81(2/ 2 +2/4)-642/ 3 ]± ¥ 8 T -R 6 . . . . 

A=252()[ 151 (2/o+2/i 2 ) +2, 048(2/ 3 +2/ 9 ) +2,4962/ 6 -2,187(2/4+2/ 8 ) 

±.022,9i2 8 

OC 

A=2^[53(y i> +y, !! )+256(y 1 +y a +y !i +yr 1 +y a +y u )+81(y<,+yi ) 
+782/ (! +109( 2 / 4 +y 8 )]±.000,14i2 8 



— 2(2/4 +2/ 8 )]±.000,099i2 8 

A =-|q[ (2/0 +2/2+2/4 +2/ 8 ) +5(2/, +2/5) +62/3] ±.006,62^ (Weddle's 
formula) 

0C 

A=^2 [17(2/o +2/t 8 ) +189(2/ 2 +2/, 0) +297(2/4 +2/ 8 )-256(2/ 3 +2/9) -742/ 6 ] 

±.002,55i2 8 (9r).* 

"Note that ordinates are not all equidistant, but would be so if certain missing ordinates. as indicated by the 
subscripts, be supplied. 



(52). 

(53). 
(53a). 

(54). 
(54a). 

(55). 
(55a). 

(56). 



(lr). 
(2r). 

(3r). 
(4r).* 

(5r).* 
(6r). 
(7r). 
(8r). 
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A=gJ^[217(2/ +2/ 8 ) +1,024(2/, +y s +2/5+2/7) +S52(y, +y 6 ) +436|/ 4 ] 

±.002,02i2 8 (14r). 

A~ [41(y +y (i )+216(y 1 +y 5 )+ZI(y,+y i )+Z12y i ]±.008,39R 8 (12r). 

A^ Q [15(y +y 6 )+U(y 1 +y 3 +y 5 )+29(y i +yJ]±.006,05R e . (18r). 

A=^ l25(y + y 9 )+Sl(y 1 +y 2 +y l +y 5 )+46y 3 ] ±.017,1# 6 . . (19r). 

A=2^o [49(yo +2/i 2) +288(2/! +2/ 5 +2/7+2/1 1) +448(i/ 3 +2/ 9 ) -27(i/ 2 +2/, „) 

-63(2/4+2/ 8 ) +1342/ 6 ]±.000,071 ( 2E lo . . . . . (17r). 

A- I gQ^ ) [41,833(2/o +2/, ,) +248,832(2/! +2/5 +2/7 +2/i 1) +395,264(2/3 +2/ 9 ) 
+118,4162/ 6 -29,160(2/ 2 +2/ 10 )-63,909(2/4+2/ 8 )] ±.000,095,622,2 (16r). 

Formulae expressed in terms of ordinates and the derivatives 
of the extreme ordinates, wherein E n = \_f n {a+x) —/"(a)]. 

A=Jq[7(2/o+2/*)+162/J--|^ ±.115726 (25s). 

A-~[13(|/ +2/a)+27(2/+2/ 2 )] -^^±.019,1^ . . . (26s). 

A=j^[217 (2/0+2/4) +512(2/+2/3)+4322/ 8 ]-J- 2 E,±.005,9i2 8 . (27s). 

A= I g- [239(2/ o +2/e)+729(2/ 2 +2/4) -2562/ 8 ] -j^,±.019 t 6JB 8 . (28s).* 

A^~^[7, 435(2/o +2/, .) +32,768(2/3 +2/ 9 ) +34,560i/ 6 -19,683(2/4+2/ 8 )] 

-J-^±.002,88E 10 (29s).* 

A^jgl^^^d/o ±i/ 12 ) +3, 402(2/ 3 +2/, 0) -2,048(2/3 +2/ 9 ) +3,888(2/4+2/8) 

+l,3402/ 6 ] -^ 2 ±.000,151i2 10 (32s).* 

A =j^ [3,149(i/„ +2/ 6 ) +7,776(2/, +2/5) +6, 075(2/ a +2/4) +8,0002/.,] 

-^J? 2 ±.000,633-R, o (35s). 



JL 
600' 

x 



A=^gJgQQQ- [13,929,421(2/0+2/, 8 ) +35,831,808(2/ 1 +2/ 5 +2/7+2/1 1 ) 

+25,824,096(2/ 2 +2/io)+38,174,720(2/ 3 +2/ 9 ) +25,511,355(2/4+2/ 8 ) 

•Note that ordinates are not all equidistant, but would be so if certain missing ordinates, as indicated by the 
subscripts, be supplied. 
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+28,171,58% 6 ]-jJ- # 2 ±. 000, 002,872, 4 .... (39s). 

A= I | 6 [31(y +2/ 2 )+6% 1 ] ~^ 2 E,+ —^E^.mARs . . (42t). 

A ~ g [121 (i/„ +2/3) +243(2/, +y 2 )] -|jg #.+ g|^#4±.006,03i2 8 (43t). 

A=3^30[3,937(2/o+2/4)+8,192(2/ 1 +2/ 8 )+7,8722/ 2 ] -^ #* + 25?040^ 4 

±.001,3172, (44t). 

A= 1635480 £ 132 ' 761( ^° + ^ s) +279,936(2/ 1 +2/ 6 )+264,627(y 2 +^) 

+280,832^] - -g^- ff a + 1401840 ff4 ± -000,07972, , . . (50t). 



The well known Simpson's rule and Weddle's rule involve the succes- 
sive use of Simpson's formula and Weddle's formula respectively. (See 
formulae (lr), (2r) and (8r) above). Similar use may, of course, be made 
of other of the foregoing formulae in y. 

Still other formulae may be had in these ordi nates by eliminating be- 
tween two or more equations involving the same set, as between (12r) and 
(18r). The new formulae thus formed are usually less accurate than those 
above given, having been formed by considering the remainders zero or 
negligible, but may be advantageously used in cases where the values of 
certain ordinates are not readily obtainable. 

For ready reference the numerical coefficients in some of the forego- 
ing formulae and their logarithms are given below in decimal form. There 
may be an error of ± 1 in the last decimal place, as results were obtained 
by using a 10-place table of logarithms. 
Formula (15) 



C s 


3.842,785,222 


.584,646,111,9 


c e 


3.570,612,245 


.552,742,690,0 


c. 


.866,878,308 


9.937,958,135,2—10 


C 3 


.142,012,987 


9.152,328,062,0-10 


c t 


.002,962,963 


7.471,726,222,9—10 


c, 


.000,001,260 


4.100,288,905,5-10 


Formula (21) 






^] 6 


1,450,463,049 


.161,506,669,3 


C 8 


.481,599,060 


9.682,685,630,1-10 


c 4 


.031,604,938 


8.499,754,946,5—10 


c. 


.000,470,312 


6.672,385,673,4-10 


c, 


.000,001,383 


4.140,906,756,4—10 
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Formula (22) 




C24 


1.686,455,775 


C12 


.994,333,666 


Cs 


.384,278,522 


c e 


.079,346,939 


c 4 


.003,095,994 


c s 


.000,150,278 


c 2 


.000,000,592 


c, 


.000,000,000 


Formula (24) 




O1 8 


1.814,431,370 


c 9 


.832,473,483 


c 6 


.018,318,060 


c 3 


.000,286,220 


c a 


.000,010,279 


c, 


.000,000,005 


Formula (39) 




Cl2,l 


1.136,381,334 


^6,1 


.158,693,878 


C44 


.024,767,952 


Cz,l 


.002,479,592 


C"2,l 


.000,024,187 


Ci,i 


.000,000,003 


Formula (50) 




C6.2 


1.026,986,573 


^3,2 


.028,081,664 


C2,2 


.001,095,703 


Cl,2 


.000,000,611 


Formula (51) 




Cs,2 


1.020,853,761 


C42 


.020,935,477 


C2,2 


.000,081,779 


Cl,2 


.000,000,061 


Formula (52) 




Cl2,2 


1.033,726,545 


C$,2 


.036,089,574 


Ci,i 


.002,503,394 


C3.2 


.000,140,975 


Ci,2 


.000,000,611 


Cl,2 


.000,000,000 



.226,974,956,5 
9.997,532,144,1-10 
9.584,646,111,9—10 
8.899,530,176,2-10 
7.490,800,103,7-10 
6.176,896,253,7—10 
3.772,322,141,0-10 
9.185,285,023,3-20 



.258,740,545,7 
9.920,370,408,9-10 
8.262,879,471,3-10 
6.456,699,497,3-10 
5.011,971,771,6—10 
1.737,981,960,5-10 



.055,524,091,1 
9.200,560,171,9-10 
8.393,890,090,8-10 
7.394,380,197,9—10 
5.383,590,134,1-10 
1.422,074,122,7-10 



.011,564,765,5 
8.448,422,840,3-10 
7.039,692,767,8-10 
3.786,354,762,5-10 



.008,963,532,4 
8.320,882,870,6—10 
5.912,642,905,3-10 
2.784,243,636,6-10 



.014,405,668,3 
8.557,381,757,7-10 
7.398,529,158,3-10 
6.149,141,792,4-10 
3.786,169,210,6-10 
9.222,593,207,8-20 
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Examples. 

for all in- 

X 

tegral values of n from 1 to 10, inclusive. 

Solution. I dx= 2 I — - — : — '-ax. Similarly for n=9, 

8, 7, etc. 

Each value to be found may therefore be derived by the aid of the 
sum of all the preceding values, already found. 

To determine which of the given formulae is sufficiently accurate and 
at the same time involves the minimum amount of labor is now advisable. 
Formula (17) has ±.0000712S IO for its remainder term, and, as seems prob- 
able, ought to be within the bounds of accuracy required. Expressed more 
fully this remainder term is 

712 f ' ° — — — — 

± TK 7 .j^[f' (x^)—f' (x s )], where <», < * and < x s < *. 

Developing sinx/a; into a power series in x, and differentiating nine 
times to obtain the tenth derivative, we get another power series in x, the 
value of which lies somewhere between and — | for < x < *. Substituting 
- ^ for the bracket expression above, we find that the remainder term has a 
value within the limits of ±.000,000,4. (17) is therefore sufficiently accur- 
ate for this value of x, and it can be readily shown that its value for other 
values of x involved in the problem is less than its value between the limits 
here considered. 

The chordal areas (represented generally by C with a subscript equal 
to the number of equidistant ordinates employed, less one, the first ordinate 
being taken at the point x=0 and the last at the point x=X, where X is the 
maximum value of x), involved in (17) are, 

C, 4=1.850,119, C e =1.844,651, C 4 =1.835,508, 
C s =1.822,636, and C 2 =1.785,398. 

Substituting in (17), J" —dx =1.851, 937. 



To check the accuracy of the numerical work done, these values of C 
may be substituted in (8) , and, by doing so, it is found that the same result 
is obtained by either formula. It is therefore probably safe to employ the 
simpler formula in computing the remaining terms sought. 

It is also instructive to apply some of the other formulae involving 
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the same values of C, e. g., (9), (10), (11), (13), (16), etc., as well as 
those involving the second, fourth and sixth derivatives of the function 
sought, such as (43a), (45a), (48a), (55a), etc. 

Ex. 2. Find the common logarithm of 2 from the identity, 



101og2-logl000=ilff 



• 024 dx 



o 1+x' 

M being the modulus of the common system of logarithms. 

Here the remainder term is easily evaluated, and the selection of the 
formula to be used is determined largely by the degree of accuracy required. 
Employing (35s) we get, log2=.301,029,995,663,981,195,213, 738,269,3, which 
we know from the form of the remainder must be correct to not less than 21 
decimal places. Expressed in terms of the formula after some simple reduc- 
tions have been made, the development is, 

l q^ * , M/ 253 [" 3,149 7,776 . 6,075] 8 2,277 
log z to + 1Q ^500 [HMO 64,005 + 64,008 J 1,771 10 8 .512 

±.000,633i? 10 . 

Using (39s) a result correct to thirty places of decimals would have 
been assured. Should still greater accuracy be desired, the step-by-step 
process of Ex. 1 may be used. 

Ex. 3. Similarly, the logarithm of 2 known, the common logarithm Of 
3, correct to forty-seven decimal places, may be obtained by (39s) from the 
identity, 

81og3+loglO-161og2=M f " Vr?-. 

Ex. 4. Obtain the value of ( - — ... , , rrr correct to six deci- 

J o sinx(l+Acos 3 a;) 2 

mal places. 



